Spectroscopy of Superfiuid Pairing in Atomic Fermi Gases 
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We study the dynamic structure factor for density and spin within the crossover from BCS super- 
fluidity of atomic fermions to the Bose-Einstein condensation of molecules. Both structure factors 
are experimentally accessible via Bragg spectroscopy, and allow for the identification of the pairing 
mechanism: the spin structure factor allows for the determination of the two particle gap, while the 
collective sound mode in the density structure reveals the superfiuid state. 



Atomic fermions have attracted a lot of interest as cur- 
rent cooling techniques allow for the creation of molecular 
condensates |lj, |2j, Li L| • These superfluids behave very 
much like standard Bose-Einstein condensates (BEC): 
the condensate may be inferred from the momentum dis- 
tribution measured in a time of flight experiment. The 
tunability of the interaction through a Feshbach reso- 
nance then offers the possibility to explore the crossover 
from BEC of tightly bound molecules to the BCS super- 
fluid state, where Cooper pairs only exist due to many 
body effects 0, 0, R, |M Ifll JlfTl. Recent experiments have 
entered this regime |lllll2Ul3j ]. however, clear signatures 
for extended Cooper pairs in a BCS like ground state are 
missing so far. In this letter, we present a generalization 
of available spectroscopic tools to measure the dynamic 
structure factor for density and spin, which reveal impor- 
tant information on the pairing mechanism within the 
BEC-BCS crossover. 

In conventional superconductors, the main character- 
istic properties are dissipation free transport and the 
Meissner-Ochsenfeld effect, which reveal themselves in 
the current response; density fluctuations are suppressed 
due to long-range Coulomb interactions 0] (see Ref. § 
for the current response in the BEC-BCS crossover). In 
contrast, for uncharged atomic gases transport measure- 
ments are not readily accessible due to the trapping po- 
tential. Then, the dynamic structure factors for den- 
sity and spin are suitable quantities for the characteriza- 
tion of the superfiuid ground state within the BEC-BCS 
crossover. Both quantities are accessible in traps: recent 
experiments measured the dynamic structure factor in in- 
teracting Bose gases via Bragg spectroscopy 0, 0, ^| , 
while the dynamic spin susceptibility may be inferred by 
measuring the spin flip rate in stimulated Raman tran- 
sitions (la . Il9j |. In this paper, we analyze the dynamic 
structure factor S c and the dynamic spin structure fac- 
tor S s within the BEC-BCS crossover. We find, that the 
dynamic spin structure factor is dominated by processes 
which break paired fermions into two single particles and 
therefore reveals the many-body excitation gap. Further- 
more, it provides the density of states, which signals the 
BCS pairing mechanism via the the appearance of a van 
Hove singularity. The observation of this singularity was 
a fundamental indication for the validity of BCS theory 



in conventional superconductors [14]. In turn, the su- 
perfiuid transition is characterized by the appearance of 
a collective sound mode in the dynamic structure factor; 
this collective mode has recently been studied [IflEfll2ll |. 

An interacting atomic gas of fermions with density 
n F = fcp/37r 2 and two different spin states is charac- 
terized by the scattering length a F allowing to tune 
the BEC-BCS crossover via the dimensionless parameter 
l/(fc F a F ). As shown by Nozieres and Schmitt-Rink [jj, 
the BCS wave function becomes exact in the BCS limit 
(l/(fc F a F ) < -1) and the BEC limit (l/(fc F a F ) > 1). 
In the following, we use this BCS wave function to de- 
termine the dynamic structure factor S c from the den- 
sity response function % c , and the dynamic spin struc- 
ture factor S s from the spin susceptibility Xs ■ This wave 
function accounts for the pairing between two fermions, 
while the residual interaction between unbound fermions 
providing particle-hole scattering is neglected. The in- 
teraction between molecules (Cooper pairs) is accounted 
for within Born approximation @ (it is not necessary to 
introduce a molecular field as done in Ref. [lfj). Within 
the BCS variational wave function, the fermionic normal 
Green's function Q and the anomalous Green's function 
T + = —T take the standard form 01 
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Here, E^ = \J t\. + A 2 denotes the single-particle excita- 
tion energy with ek = fi 2 k 2 /2m-/i the free fermionic dis- 
persion relation, while O s = ttT(2s + 1) denotes fermionic 
Matsubara frequencies. The presence of a condensate 
and superfiuid response in the ground state is encoded 
in a finite BCS gap A. The gap A is determined by 
the scattering length a F via the gap equation, and the 
chemical potential is fixed by the particle density n F 0, 
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In the BCS limit at low temperatures (T = 0), these re- 
lations give the chemical potential fi — e F and the gap 
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A = e F (8/e 2 ) exp (7r/2fc F a F ) with e F the Fermi energy. 
In turn, in the BEC limit the appearance of a two par- 
ticle bound state with binding energy e = SVW) 
and internal wave function 0k = (2n F ) -1 / 2 A/E^ mod- 
ifies the chemical potential p = —cq/2 and the gap 



Within the crossover regime, we solve 



Eq. and (0) numerically to determine A(fc F a F ) and 
p(k F a F ). Note, that in general the BCS gap A differs 
from the two particle excitation gap A s = min k 2Ek. 
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FIG. 1: Diagrams contributing to the density response func- 
tion and the spin susceptibility; Xpair accounts for particle- 
hole excitations, while the diagrams in Ycoii describe collective 
sound excitations. 

In the following, we calculate the dynamic structure 
factor S c — — Imx c /7r and the dynamic spin structure 
factor S s = — Imxs/7r via their relations to the density re- 
sponse function and the spin susceptibility, respectively. 
The density response 5p c {uj, k) to a small external drive 
6V c (u!,k) follows from <5p c (a;,k) = Xc(<^, q)SV c (uJ, k) 
with the linear response function (in real space) 



xo(*,x) = -*e(t)<|> (t,x),p (o,o)]>. 



(3) 



Here, (. . .) denotes the quantum statistical average at 
fixed temperature T and chemical potential fx, while 
the density operator is defined by p c = ip^ip^ + V'J'V'j.- 
The analogous definition applies for the spin suscepti- 
bility with the density pc replaced by the spin density 
p s = if>^if)f — The diagrams contributing to the 

response functions Xc and Xs are shown in Fig.^ Note, 
that each diagram has to be weighted by a nontrivial 
factor ±1. These factors differ for the response func- 
tion Xc and the spin susceptibility Xs, and provide dif- 
ferent cancellations between diagrams. We distinguish 
between two types of diagrams: The first type of dia- 
grams Xpair involve only normal and anomalous Green's 
function and describe two particle excitations, while the 
second type of diagrams x&>n involve the vertex opera- 
tor T a /3 and account for collective excitations. The BCS 
wave function neglects particle-hole scattering and T ac- 
counts only for particle-particle (hole-hole) scattering. 
Then, a,/3 <E {p,p} describe the incoming and outgoing 
type of particles; p accounts for particles and p for holes. 
The vertex operator T has to be calculated with the help 
of the BCS wave function, see below. Note, that Fig. 
shows only the diagram involving T pp , while the other di- 
agrams exhibit a similar structure. In the following, we 
are mainly interested in the zero temperature limit T = 



and in the low momentum regime k <C l/£ = A/hc s (the 
energy of the collective modes is below the two particle 
gap uj < A s ). Here, c s denotes the macroscopic sound 
velocity and £ the size of the pairs with £ ~ hv F /A in 
the BCS limit and £ ~ a F in the BEC limit. For typi- 
cal experimental setups, the scale £ is small compared to 
the trap size R, and we can safely assume the condition 
1/R < k < l/£. Then, the trapping potential plays a 
minor role as has been shown in the measurement of the 
dynamic structure factor and the tunnelling probability, 
see Ref. 0] . Therefore, we ignore the influence 

of a trapping potential in the following analysis. 

First, we focus on the spin susceptibility Xs- At zero 
momentum, the spin susceptibility Xs( w jO) is equiva- 
lent to the response driven by the spin flip Hamiltonian 
H = X(t) J dx\ipfip, + xj}fxj}A. A perturbation of this 
form is realized experimentally by driving a stimulated 
Raman transition between the two hyperfine states of the 
two component Fermi gas |lPJ | . Within the diagrammatic 
expansion of Xs, the diagrams in Xoii cancel each other 
and only pair excitations x P air survive; their contribution 
takes the form 



3 (tt s ,k)=-2T^ 
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g(n t ,q)g(n s+t ,q + k) 



-T+(n t ,q)T(n s+t ,q + k) (4) 



(Note, that for the density response function Xc the + 
sign between the two terms is replaced by a — sign (l^.l 
The integration in 10} involves standard methods, and we 
present here only the final result for the spin structure 
factor in the low momentum limit k <^ l/£, see Fig. El 
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The spin structure factor exhibits a gap A s = 2 min k E^, 
i.e., for positive p, the spin gap A s = 2 A is determined 
by the BCS gap, while for p < 0, it is given by the bind- 
ing energy A s = 1\J p 2 + A 2 approaching e in the BEC 
limit. The shape of the spin structure factor differs in 
the two limiting regimes: the shape exhibits the char- 
acteristic 1 / \Ju> — 2A singularity of the density of states 
for the BCS pairing mechanism, while in the BEC regime 
S s exhibits a maximum. In the crossover regime, Eq. JSJ) 
smoothly interpolates between these two limits. Next, we 
analyze the modifications of the spin structure factor for 
temperatures above the superfluid transition tempera- 
ture T c . Then, Eq. (JIJ implies A = 0, and the spin struc- 
ture factor in the BCS limit reduces to that of a Fermi 
gas. In turn, the presence of bound fermion pairs dom- 
inate the spin structure factor even above T c (T < T*) 
and S B exhibits a pseudo gap ~ A s eo with disso- 
ciated fermion pairs providing a small but finite weight 
within the gap. The pseudo gap finally disappears above 
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the pairing temperature T > T* via a smooth crossover. 
Therefore, the measurement of the spin structure factor 
provides a suitable tool for the characterization of the 
pairing temperature T* and the spin gap A s . 
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FIG. 2: Spin structure factor Ss(a>) in units hn F /E F for dif- 
ferent scattering lengths l/(opfep) = —0.5,-0.25,0,0.25,0.5. 
The spin structure is quenched below the spin gap uj < A s , 
and exhibits a characteristic singularity in the BCS limit. 

In contrast, the dynamic density structure factor S c in 
superfluids is dominated by a collective sound excitation 
representing the Goldstone mode of the broken symmetry 
(in conventional superconductors, Coulomb interactions 
lift this mode to the plasma frequency). The dynamic 
structure factor at low momenta k <C l/£ takes the form 

hk 

S c (w,k) = n F - 5(u — c s k) (6) 

2mc s 

exhausting the /-sum rule and compressibility sum rule, 

h 2 k 2 
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The determination of the sound velocity within the BEC- 
BCS crossover requires the calculation of the diagrams 
Xcou m Fig- Ql its contribution exhausts the sum rules 
Q and © at small momenta k <C l/£ and frequencies 
hio < A s i.e., a cancellation appears between x P air and 
Xcoii at frequencies uj > A s . The diagrams in Xom provide 
the response function 

Xo(w, k) = 4^ ^(w, k)r a>/3 (w, k)^(w, k), (9) 

where one factor 2 accounts for summation of spin in- 
dices, while the second factor 2 appears as each vertex 
operator T a p contributes to two diagrams. Using the 
microscopic approach by Haussmann |jj, the vertex op- 
erator T a /3 takes the form 



with 



1 _ 


m 


f dk m 




Anh 2 a F j 


' (2Tr) 3 h 2 k 2 



M, 



i>i> 



M, 



i>i> 



S af 3 + M af3 (10) 

[a(n._ ts q+k)e(fi t ,-q)], 
[^(n a _ t ,q+k)^(nt,-q)], 



with M pp = Af^ and M pp - = M£ p . The propagators 0* 
(0 P ) account for the creation (destruction) of a fermion 
pair from the condensate, and take the form 



b P = T Hj (03 G(Vt, q)-F(fii + .,q + k). (11 



Solving the above equations in the limit of small fre- 
quency uj <C A s and momenta k <C l/£ provides the 
collective sound mode with sound velocity c s ; the result 
is shown in Fig. [3] 
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FIG. 3: Left: Sound velocity c s in the BEC-BCS crossover 
in units v F / \/3- The dashed line is the Bogoliubov sound ve- 
locity. Right: Structure factor Sc(u) of weakly interacting 
fermions and bosons above the superfluid transition tempera- 
ture T = e F /2 > T c and k — v^^f/40. The gray line denotes 
the peak of the sound mode below the superfluid transition. 

Within the BCS limit, the collective excitations © 
provide the Bogoliubov- Anderson sound mode for a neu- 
tral superconductor with c s = v F /^/3 22]. The particle- 
hole contributions, dominating the structure factor for a 
weakly interacting Fermi gas, are quenched due to the 
opening of the excitation gap. Note, that the leading 
correction to the sound velocity c s = y F l y/3[l — 8k F a F /ir] 
derives from particle-hole scattering |22j ; such scattering 
processes are not contained in our approach. In turn, 
within the BEC limit Eq. © gives the structure factor 
S% EC (u,k) = 2n F S(uj - Mc 2 /4m). The variational BCS 
wave function approach provides the zeroth order result 
describing a noninteracting Bose gas of molecules with 
density n B — n F /2 and mass m B = m/2. The structure 
factor becomes 4 times the structure factor of an ideal 
Bose gas. This factor 4 appears as the external poten- 
tial drives the fcrmionic density operator instead of the 
bosonic density operator; the structure factor satisfies the 
/-sum rule for fermions with n F /m = 4n B /m B . Going be- 
yond leading order, the above equations also incorporate 
the repulsion between the bound pairs, and provide the 
structure factor S° EC (u), k) = 2n F ek/i?k£ (u— Ek/h) with 
e = fi, 2 k 2 /2m B and the Bogoliubov excitation spectrum 
E^ = e 2 . + 2/i B e r k- The structure factor describes a weakly 
interacting Bose gas with sound velocity c a = ^/ [i B jva B . 
Here, /i B = 4:irh 2 a B n B /m B denotes the bosonic chemical 
potential accounting for the scattering length a B = 2a F 
within Born approximation; its exact value a B w 0.6 has 
recently been derived [23j|. 
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Next, we focus on the dynamic structure factor above 
the superfluid transition temperature T > T c and com- 
pare it with the structure factor Eq. ©. We first focus 
on the collisionless regime lot > 1 with r the collision 
time; this limit is naturally achieved for weakly interact- 
ing atomic gases (|fc F a F |< 1) and frequencies to above 
the trapping frequency |2^. Within the BCS limit, the 
system reduces to a Fermi liquid with a weak attractive 
interaction. The structure factor exhibits the particle- 
hole excitation spectrum of a weakly interacting Fermi 
gas at finite temperature (interactions only renormalize 
the Fermi velocity u F ), while the zero sound mode is 
overdamped for attractive fermions; the structure factor 
is shown in Fig. In turn, in the BEC limit the sys- 
tem above the critical temperature T c reduces to a gas of 
bosonic molecules. The structure factor of a degenerate 
Bose gas with temperatures above the superfluid transi- 
tion temperature derives from the bosonic Lindhard func- 
tion; the structure at low momenta is shown in Fig. 
Comparing these structure factors with Eq. ©, we find 
that the superfluid state is characterized by the appear- 
ance of a collective sound mode. However, for the system 
in the hydrodynamic regime lot < 1, the structure fac- 
tor is exhausted by the hydrodynamic sound mode (first 
sound) even above the superfluid transition. Therefore, 
the identification of the superfluid transition from the 
density response requires to be in the collisionless regime, 
which is reached for sufficiently weak interactions. 

Within recent experiments 0, El > the molecular bind- 
ing energy in the BEC limit was determined by a r.f. 
pulse breaking the molecules and exciting a fermion into 
a different hyperfine state. This method has the disad- 
vantage that the interactions between the fermions in dif- 
ferent hyperfine states produces non-trivial energy shifts. 
In turn, the structure factors presented here, produce ex- 
citations within the same hyperfine states. However, the 
two hyperfine states are separated by an energy fko^i 
due to Zeeman splitting. Using a driving field at fre- 
quency Lo-^i with an additional superimposed modulation 
frequency lo, i.e., A(t) ~ cos(wj|i) cos(wi), then probes 
the structure factors S s (lo) with the number of parti- 
cles in each hyperfine state conserved. Such a procedure 
avoids non-trivial energy shifts induced by a change in 
particle number or excitation of particles into different 
hyperfine state, and therefore represents a suitable setup 
for a determination of the two particle excitation gap A s . 
The measurement of the structure factors can be achieved 
in two different ways: First, the energy transfer W to 
the system satisfies W = ^loS{lo, k) with 7 determined 
by the driving field alone and allows for the determina- 
tion of S from the heating of the system 0, 0] . For 
6 Li, the characteristic parameters at the Feshbach reso- 
nance are given by A s ~ E F ~ 1/iK and wj^ ~ 80MHz, 
i.e., the structure factor is accessbile via Raman transi- 
tions or r.f. pulses as and shown in recent experiments 
HI EU E3 ■ Second, the interaction between the driving 



field and the fermions leads to the absorbtion and emis- 
sion of photons with a rate determined by the structure 
factor S(u), k). Therefore, an analysis of the counting 
statistic of the probing laser field allows for an in-situ 
measurement of the structure factors. 
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